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Abstract. The spectrum of Harper's equation is determined by the dis- 
criminant, which is a certain polynomial of degree Q if the commensurability 
parameter of Harper's equation is P/Q, where P, Q are coprime positive in- 
tegers. A simple expression is indicated for the derivative of the discriminant 
at zero energy for odd Q. Three dominant terms of the asymptotics of this 
derivative are calculated for the case of an arbitrary P as Q increases. The 
result gives a lower bound on the width of the centermost band of Harper's 
equation and shows the effects of band clustering. 

It is noticed that the Hausdorff dimension of the spectrum is zero for the 
case P — 1, Q infinitely large. 
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1 Introduction 



In the present paper we consider the following eigenvalue equation in l 2 (Z) 
(see fill HIl and [El-Sl for the reviews): 



ip n - 1 +2cos(27ran+9)ip n +ip n+ i — £ip n , n = . . . , — 1, 0, 1, . . . , a,9e~R 

(1) 

It is called Harper's equation with commensurability parameter a. If a = 
P/Q, where P and Q are relatively prime integers, the spectrum of the 
corresponding operator consists of Q intervals (bands). The interesting case 
is the incommensurate one, i.e. when a is irrational. It is known that for 
irrational a (but not for rational) the spectrum is independent of 9. There 
is a conjecture that for any irrational a the measure of the spectrum is zero. 

Let us call a(a) the union over all real 9 of the spectra of (U). Given 
(|ip with a = P/Q, one can construct an independent of 9 polynomial S(x) 
of degree Q with the following property: the spectrum a (a) is is the image 
of the interval [—4,4] under the inverse of the transform T,(x) = A. For 
convenience, we call the discriminant of Harper's equation.^] It is known 
that S(x) = — £(— x) (hence S(0) = 0) for odd Q and S(x) = £(— x) for 
even Q. Also for even Q we have |S(0)| = 4 (which implies that 2 bands of 
the spectrum merge in 0). In the present paper, we will give the expressions 
for the derivative at zero S'(0) for odd Q. Such expressions provide the 
lower bound on the width of the centermost band. Using our approach, one 
can also calculate higher derivatives at zero, but the formulas become more 
cumbersome. 

Several important results about Harper's equation with irrational a have 
been obtained using the properties of the discriminant. Still not much is 
known about these properties. Here we briefly mention some results. In 
1992, Last and Wilkinson || noted the following fact about the derivatives 
of at its zero crossings: 

Q ! ! 

where x/-, k — 1, . . . , Q are the zeros of E(x). This property was subsequently 
used by Last [IU| to prove the vanishing of the measure of the spectrum for 



1 The usual definition of the discriminant of a periodic matrix corresponding to ([!]) 
would be S(x,9) = £(x) - 2 cos 9Q. 



2 



a certain class of irrational en's. More detailed information on can be 
obtained in the semiclassical case, that is when a is small [11,1^1 (also when 
a is close to a rational [|13|| ). The uniform asymptotics is available for the 
case when P = 1 and Q is large 



141. From it, one can obtain the widths 



and distribution of individual bands. The large Q asymptotics can also be 
calculated for P = 2, 3... [fbl but the expressions quickly become cumbersome 



with growing P. As we shall see, this is not the case with the asymptotics of 
S'(0). We will calculate it for an arbitrary fixed P up to the terms decreasing 
with Q. We will obtain several different types of asymptotic behaviour, which 
reflects the effects of band clustering. 

There is a considerable amount of literature on the question of what is 
the Hausdorff dimension of a(a) for irrational a and references therein). 
We make a note on this by showing that the Hausdorff dimension for the case 
of P = 1 and infinite Q is zero. 



2 General formula for the derivative of the 
discriminant at zero 

The main fact which will enable us to obtain an expression for S'(0) is that 
the spectral problem (1) can be reduced to that for a matrix with zero main 
diagonal. In various avatars, this fact was noted in JT6| , pT| , P^[ . 



Henceforth, we assume a = P/Q, where Q is odd, P is relatively prime 
with Q. It follows, e.g., from the reasoning in Appendix of pOf that = 



— det(L — xl), where L, I are Q x Q matrices, / is the identity matrix, and 
the matrix elements of L are the following: L\. u-i = L^-x k = 2 sm(irPk/Q), 
if k = 1, 2, . . . , Q — 1, and zero otherwise. Thus, L is a tridiagonal matrix 
with zero main diagonal. It is therefore easy to get a simple expression for 
the derivative of the determinant det(L — xl) at x — 0. Using the identity 
ni=7 1 ^ 2 (2 sm(irP2k/Q)) 2 = Q at the last stage of the calculation, we obtain 

E «^-.)'-'''«(HTnf^), <3) 

where 7 = nP/Q. In the next section, we shall calculate the large Q asymp- 
totics of S'(0) for the case of fixed P. 
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3 Asymptotics for E'(0) 

Let P be arbitrary fixed and consider the asymptotics of S'(0) for large Q. We 
can represent any odd Q (relatively prime with P) in the form Q = 4Pr + s, 
where r,s are nonnegative integers, s = 1,3,5, ...,4P — 1, excluding the 
numbers which are not relatively prime with P. Let 

(Q-i)/2 k s i n 2(2 ? -lU 

S= V TT [ 4 (4) 

tx i=i ^n 2 2j 7 

and rewrite it in the form 

(Q- ! )/4P 2 p-i^+- j 2f2 _ ^ ( s -i)/2 sin 2 (2 7 - 1W 

s= E I] TT 1 o n + E TT 1 o n (5) 

^ h tA sin 2 2 J7 ^ ^ sin 2 2 J7 1 ' 

To obtain the asymptotics, we first break the products here into subproducts 
with the number of factors (Q — s)/4P. Then the sinuses become small for 
large Q only in one connected part of each of the subproducts. By a change 
of product variable j = k(Q — s)/4P + f for k even, and j = (k + 1)(Q — 
s)/4P — j' + 1 for k odd, we obtain: 

^ +m sin 2 (2j-l) 7 ^ sin 2 (2j- 1 + ^)7 fR . 
11 ^TTT- = A A 1 ■ ■ ■ A t ^ [[ — . 2 -— r 6 

sin 2j7 sm (2j + 5 2 )7 

where ^ = 5 2 = -ts/2P for t even; 6t = (t + l)s/2P, 5 2 = -1 + (t + l)s/2P 
for t odd; and 



(Q Y| /4P sin 2 (2j - 1 - ks/Ph 
2k ~ }} x sin 2 (2j- ks/Ph ' 

{Q-s)/AP . 2 



A 2 k-l - II , 2 



sin 2 (2j - 1 + ks/P)^ 
sin 2 (2j -2 + fcs/^b" 



(7) 



Substituting these expressions into (0) and making a change of variables 
to = (Q — s)/4P — to' in one of the sums, we finally obtain the following 
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expression that will be a starting point for deriving the asymptotics: 



5 "S5 A i £ sm 2 ( 2j -t S /P) 7 ' 

n sin2(2j _ i + (t + i )s/P ) 7 j + A ^ A ^ - x ) j + ( 8 ) 

n 1 4 ( v /2 fr sin2{2j ~ 1 ~ 5)7 

M 1 M sin 2 ( 2j - S ) 7 • 

Henceforth, we assume n°=i /» = El=i 3™ = 0. 

We will need to estimate asymptotics of 2 types of products. Let Cj = 
sin 2 (2j + 2a) 7 / sin 2 (2j + 26) 7 , where a and b are any constants such that Cj 
(j = 1, 2, . . .) are finite. First, for k(Q) such that k/Q — > as Q — > oo, we 
have: 

17 c - r 2 (& + i)r 2 (fe + a + i) oro _ 2+ ^ f9 . 

/i Cj "-p(a+l)P(]fe + 6+l) (1 + 0(g )} ' (9) 

where we can choose e arbitrary small if k/Q — > fast enough. Secondly, if 
M — > oo, but M/Q — > as Q — > oo, we have 

n * = 

j=M+l 

j * l-(2a 7 ) 2 + 0( 7 4 ) + (2a 7 + 0( 7 3 ))cot2 J7 
exp<2 ln ~ 



(26 7 ) 2 + 0( 7 4 ) + (26 7 + 0( 7 3 )) cot 2j 7 
Sln( f M 1)7 ) ^ ^ ( X + 27 (° 2 " [cot(2A; + 1)7 " cot(2M + 1)7 1+ 



j=M+l 

2(o-6) 



2 7 M 

0(1/Q 2 ) + 0(1/MQ)). 

(10) 

Here we expanded the logarithm in series using the fact that 1/M is a small 
parameter, and then applied the Poisson summation formula to estimate the 
sum. 



5 



Now let us estimate the asymptotics of A^. 

it sin 2 (2j - 1 - ks/P)^ W-*)/ 4P S in 2 (2j - 1 - A;s/P) 7 
2k ~ Al sin 2 (2j-^/P) 7 . sin 2 (2j - ^/P) 7 

T 2 (l-A;s/2P) , , „ NN 

2 ^ P(l/2-fa/2P) (1+ ° (1/g)) - 



:n) 



where we used fllQl) and (|]) for b = —ks/2P, a = b — 1/2 (in @ we took the 
asymptotics of T(x) for large x). 

Similarly, using @ and ( |10D for 6 = — 1 + ks/2P, a = b + 1/2, we get: 

^^fTO^P) (1 + 0(1/<32)) ' (12) 

The last nontrivial step in our calculations is the following estimate for 
a = 6-1/2: 



?-s)/4P fc M fc (Q-s)/4P M fc 

E II o Ello E Ho II 

fe=l j=l k=lj=l k=M+l j=l j=M+l 

f2(6 + 1} >ri(b) - In 7 + (4P - s)^- + 7 (26 - ^) + 0(1/Q 2 ^ 



(13) 



r 2 (6+ 1/2) [ ,w ' v y 2g ' V 2 

where the constant of Euler's type 

" ' r 2 (t + 6+i) / ^ ' 

Such constants were studied in 

To obtain (p~3| ) we used expressions (|9j|10l), and then applied the Poisson 
summation formula to get the asymptotics of the sum. The most important 
moment in our derivation is to choose M in such a way that M — *■ oo and 
M/Q — > as Q — > oo. Apart from this condition, M is arbitrary. Naturally, 
M cancels in the final expression. 

Using (|13|) pTl , p!2D and (fj), it is easy to get the asymptotic formula for 
S from (El). As a result, we find that the derivative of the discriminant 



2 I am grateful to R. Askey for indicating this reference. 



6 



E'(0) = (-1)(°~ 1 )/ 2 (5(1 + S) has the following asymptotics as Q -> oo: 

p-i t 



s'(o) = (-i)^ 1 )/ 2 g^Ericot 

's - 1 



7TS2 



2P 



2 In 



t=0 i=l 

st_\ T\s/2-st/(2P)) \ 
2 2Pj~ t T 2 ((s + l)/2-st/(2P))r 



Q_ 

7tP 



St 

V[ -2P 



(15) 



1 



(«-l)/2 r2 



7T 



E 

m=l 



r 2 (m- 1/2) 
T 2 (m) 



o(l) 



where the constants r)(b) are defined in Qlj 

It is known that the bands for large Q appear away from zero energy in 
clusters of P bands separated by wide gaps. In a cluster around zero energy 
there are s bands. We see from ( |I5] ) that at least to the leading Q In Q order 
there are no more than P/2 types of different asymptotic behaviour (because 
of the periodicity of cotangent). Namely, for odd P > 1, s = 1,3, . . . ,P — 2; 
for even P, s — 1, 3, . . . , P — 1 excluding the values not relatively prime with 
P. It would be interesting to prove if this is the case for the whole asymptotic 
series. 

For P = 1 there is only one type of the asymptotic behaviour: 



]'(0) = l(-l)(Q-WQ (in 9l + rj(Q) + n ) + (i) 

7T \ TT / 



(16) 



4 Hausdorff dimension 



Consider the case P = 1, Q — > oo. We are interested in the Hausdorff 
dimension of the limiting spectrum (there is no reason for it to be equal to 
the dimension for the case P/Q = 0, which is, obviously, 1). The spectrum of 
Harper's equation lies within the interval [—4,4]. The dominant term of the 
uniform asymptotics for the discriminant on x G (—4, 4) as Q — > oo can 
be found in |14|. (One can also give a simpler derivation of ( jI7| ) taking the 



representation of S(x) in terms of det(L-xI) and matching the semi classical 
solution with the asymptotics for Meixner-Pollaczek polynomials p2].) 

For \x\ < I, (A = x/A) 



4 cosh(A(y) cos In 2 arg Y - + % (17) 

\7T 7T \2 IX I 2 ) 
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For xQ>l, x < 4 - e < 4 (for any e > 0) 

E(x) ~ 2e 2Q ^ cos(2Qz/) 



where 



1/2 

(27T)- 1 arccos(2A-l) 
(27T)- 1 arccos(2A-l) 



(18) 



(19) 



arccosh(2A — cos27ct)dt, 
arccos(2A — cos2nt)dt, 



the branch of arccosine being < arccosa; < n for x £ [—1,1]. Given fllq) , 
the asymptotics for the symmetric region of negative x is obvious because of 
the property S(a?) = (— 1) Q S(— x). 

It is interesting to compare S'(0) obtained from ( |T7| ) with (|TB]) . This gives 
the representation 7/(0) = In 16 + C — 7T, where C is Euler's constant. 

In the region where \xQ\ ^> 1, the bands of the spectrum are exponentially 
narrow in Q, and their density, as one can deduce, is given by 



p{x) = K'(x/A)/(2tt 2 



(20) 



where K'(X) = K(^l^X 1 ) = / " /2 (l - (1 - A 2 ) sin 2 u)- l / 2 du is the complete 
elliptic integral. Indeed, we see from ([Tj^) that p(x) is given by the expression 
p(x) = (l/7r)2\dv/d\\d\/dx. Further, we have 



dv 




d\ 


= V 

7T JO 



arccos(2A— 1) 



dt 



Making the change of variable cost 



dv 



dX 



2/ n 



1 + A U + A 



K 



1 - (2A-cost) 2 
- (1 — X)u + A, we get 
1-A\ 1 



(21) 



7T 



-K'(X), 



(22) 



which yields (0). (One can also deduce (|20|) from the results of [0]). This 
expression for p(x) as a = P/Q — > is very natural as it is exactly the 
well-known density of states for the case a = 0. 

Because of the exponentially small widths of the bands, it is easy to show 
that the Hausdorff dimension for the region \xQ\ ^> 1 of the spectrum is 
equal to zero. (Note that for x G [4 — e, 4], < e <^ 1, the bands remain 
exponentially small. The expressions for the positions of these bands to 
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several orders in 1/Q are given in JH].) In the other region, \x\ <C 1, the 
widths of the bands are of order 1/Q In Q when \xQ\ ~ 1 as we see from 
flT7fl. More precisely, the dominant asymptotic term 8(XQ) for the width of 
the bands in the neighbourhood of XQ is 

Att 1 
5(Ag) = ^ arCSm oo^hW ^ 

The factor 4 appears here because we are looking for the width "on the scale" 
x, while X — x/A. 

If we sum the widths 5,- of all these bands, we will obtain the famous 



Thouless formula for the total bandwidth W . Indeed, if t = XQ, then 

in the interval dt there are, to the main order, - lnQ— bands. Hence, 

' ' TT ^ TT ' 

16 r°° 1 32 

W = £ft~-n/ arcsin— -rft=— /3(2), (24) 
7rQ Jo cosh t tiQ 

where f3(2) is Catalan's constant. To estimate the Hausdorff dimension, we 
consider the sum W(d) = J28f- We have 

I /± l + d poo i 

W(d) ~ -r-, r-j / arcsin d — —dt, (25) 

V ; Q d \n d - 1 Qn 2 - d Jo cosht v ; 

which tends to zero as Q —>■ oo for any d > and is infinity for d = 0. This 
implies that the Hausdorff dimension is zero. Thus, the Hausdorff dimension 
of the whole spectrum for P — 1, Q — > oo is zero. 
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